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Abstract 

The electromagnetic form factors of the nucleon, in the space-like region, are determined from three- 
point function Finite Energy QCD Sum Rules. The QCD calculation is performed to leading order in 
perturbation theory in the chiral limit, and to leading order in the non-perturbative power corrections. 
The results for the Dirac form factor, Fi(q 2 ), are in very good agreement with data for both the proton 
and the neutron, in the currently accessible experimental region of momentum transfers. This is not 
the case, though, for the Pauli form factor ^(g 2 ), which has a soft g 2 -dependence proportional to the 
quark condensate < 0|qq|0 >. 
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The electromagnetic form factors of the nucleon have been studied in perturbative QCD (PQCD), together 
with QCD sum rule estimates of the nucleon wave functions [1]. Comparison with data is difficult due to 
the extreme asymptotic nature of these theoretical results. In fact, the onset of PQCD in exclusive reac- 
tions does not appear to be as precocious as in inclusive processes. In addition, these wave functions are 
affected by some unavoidable model dependency. In any case, the Dirac form factor Fi(Q 2 ) does exhibit 
the expected leading asymptotic 1/Q 4 behaviour. However, the Pauli form factor ^(Q 2 ) turns out to be 
of higher twist, and therefore not accessible in the standard PQCD approach. At current experimental 
space-like momentum transfers, the results from the standard hard-scattering approach for Fi(Q 2 ) do not 
compare favourably with the data. On the other hand, some recent light-cone QCD sum rule determina- 
tions appear to improve the agreement with data from within a factor 5-6 to within a factor of two [2]. 
The source of this persistent disagreement does not seem easy to identify. In view of this, it is desirable 
to attempt a QCD sum rule determination in a region of experimentally accessible momentum transfers, 
and without any reference to the concept of a wave function. In addition, one should employ sum rules of 
a type which would provide a clear insight into the source(s) of potential disagreement with experiment. 
This can be achieved e.g. by using Finite Energy Sum Rules (FESR). In fact, in this framework the power 
corrections involving the vacuum condensates decouple to leading order in PQCD. In other words, power 
corrections of different dimensionality contribute to different FESR. 

In this note we determine the Dirac and Pauli electromagnetic nucleon form factors, in a wide range of 
(space-like) momentum transfers, in the framework of three-point function QCD-FESR of leading dimen- 
sionality. As is well known by now, this technique is based on the Operator Product Expansion (OPE) 
of current correlators at short distances, and on the notion of quark-hadron duality [3]. Analyticity and 
dispersion relations connect the QCD information in the OPE to hadronic parameters entering the corre- 
sponding spectral functions. We compute the QCD correlator to leading order in perturbative QCD in the 
chiral limit (m u = m d = 0), and include the leading order non-perturbative power corrections proportional 
to the quark-condensate and the four-quark condensate (with no gluon exchange). We begin by considering 
the following three-point function (see Fig. 1) 

Il^p 2 ,p' 2 ,Q 2 ) = i 2 fd 4 xfd 4 ye^' x -^ (0\T{ m (x)J* M (y)fjN(0)}\0) , (1) 

where Q 2 = —q 2 = — (p' — p) 2 > is fixed, and 

Vn(x) = e abc [u a {x){C la )u\x)] (j 5 rd c (x)) (2) 

is an interpolating current with nucleon (proton) quantum numbers; the neutron case u <-> d will be dis- 
cussed at the end. In Eq.(l), J^ M is the electromagnetic current 
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J£m(v) = \u{y)-fu{y) - \d{y)^d{y) . (3) 
The current Eq.(2) couples to a nuclcon of momentum p and polarization s according to 

(0\r lN (0)\N(p,s)}=\ N u(p,s), (4) 

where u(p, s) is the nucleon spinor, and Aat, the current-nucleon coupling, is a phenomenological parameter 
a-priori unknown. This parameter can be estimated, e.g. using QCD sum rules for a two-point function 
involving the currents ijn [3]- [4]. In this case one can determine the nucleon mass, as well as the coupling 
\n. 

Concentrating first on the hadronic sector, and inserting a one-particle nucleon state in the three-point 
function (1) brings out the nucleon form factors Fi(q 2 ), and F 2 (q 2 ), defined as 



(JV(fci , si) \JZ M (0)\ N(k 2 , s 2 )> - u N {k u si) 



u N {k 2 , s 2 ) , (5) 



where q 2 = (k 2 — ki) 2 , and k is the anomalous magnetic moment in units of nuclear magnetons (k p = 1.79 
for the proton, and K n = —1.91 for the neutron). The form factors F\ t2 (q 2 ) are related to the electric 
and magnetic (Sachs) form factors Ge{<1 2 ), and Gm(<Z 2 ), measured in elastic electron-proton scattering 
experiments, according to 

G E (q 2 ) EE Fl { q 2 ) + ^F 2 {q 2 ) : (6) 
G M (q 2 ) EE F.iq 2 ) + KF 2 (q 2 ) , (7) 

where G P E (0) = 1, G P M (0) = 1 + n p for the proton, and G£(0) = 0, G n M {0) = K n for the neutron. Next, 
the hadronic spectral function is obtained after inserting a complete set of nucleonic states in (1), and 
computing the double discontinuity in the complex p 2 = s, p' 2 = s' plane. For s, s' < 2.1 GeV 2 , i.e. below 
the Roper resonance, one can safely approximate the hadronic spectral function by the single-particle 
nucleon pole, followed by a continuum with thresholds sq and s' (sq, s' > M^). This hadronic continuum 
is expected to coincide numerically with the perturbative QCD (PQCD) spectral function (local duality). 
This procedure is standard in QCD sum rule applications, and leads to 

Imll^s, s', Q 2 ) = n 2 \ 2 N 6{s - M 2 N )8{s' - M 2 N ) 

HAD 

X \ Fl (q 2 ) [p"^ + M N {j>' ltl + + M 2 Nltl ] 
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^-F 2 (g 2 ) [fa^ + M N (fa^ + a^f) + M^a^] q v 1 6 (s - s) 
Imn M ( S ,s',Q 2 ) e(s-so) , 



PQCD 



where we have set s = s' for simplicity. 
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Turning to the QCD sector, the three-point function (1) to leading order in perturbative QCD, and in the 
chiral limit, is given by 



W(p 2 ,p' 2 ,Q 2 ) = 16 f d 4 .x f d 4 y jlp'-*-*-v) Tr [ f i*^fc e -*M*-i/ 
J J \_J (2tt) k\ 



d k 2 fa lk2 .y 



(2tt)* k\ 



d4fc 3 fa 0+ ik 3 



(2tt)4 kl 



7 5 7" 



d 4 fct fa 
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+ 4 



dye 



i(p'-x-q-y) 



Tr 



x 7 



(2tt)4 fc| 

/• d 4 fc 4 jfa-s /" d 4 fc 3 ^3 - fc3 . x 
7 (27r)*fcf 7 "y ^fcf 
d 4 fc 2 £ 



7o 



(9) 



d fci fa ikl . (x _ y) 
[2^k\ 7 



(2tt) 4 k 2 



2 e ~ik2-y^v 



After computing the traces and performing the momentum space integrations, Eq.(9) involves several 
Lorentz structures analogous to those entering the hadronic spectral function Eq. (8). Before invoking 
duality one needs to choose a particular Lorentz structure present in both (8) and (9). A convenient 
choice turns out to be ^'7^, which allows to project Fi(q 2 ), as this structure does not appear multiplying 
F2(q 2 ) in Eq.(8). An additional advantage of this choice is that the quark condensate contribution, to be 
discussed later, does not involve the structure ^'7^, on account of vanishing traces. There is, though, a 
non-perturbative term involving this structure and proportional to the four-quark condensate. However, 
eventually this term will not contribute to the FESR as its double discontinuity vanishes. Hence, Fi(q 2 ) 
will only be dual to the PQCD expression. It must be pointed out that the PQCD spectral function con- 
tains the structure f'^i^f explicitly, as well as implicitly, i.e. there are terms proportional to this structure 
which are generated only once the momentum-space integration is performed. 
After a very lengthy calculation, the imaginary part of Eq.(9) is given by 



ImIP(s, s',Q 2 ) 



(27T)' 



: (3fii + 4r! 2 -fi3) 



(M) + ■■■ 



(10) 



where 



«i = y \Q 2 + s 



Q i + 2Q 2 s + s 2 -2s s' -s' 2 
^Q A + (s- s'f + 2Q 2 (s + s')J 



(11) 
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n 2 



+ 



7T 6 < 



(2 Q 2 + 3 s - 3 s') [(Q 2 + s) 3 (2Q 2 + 3s) + 3 (Q 6 - 5QV - 4s 3 ) «'] 

Q 4 + (s-s') 2 + 2Q 2 (s + s') 



(3Q 2 - 4s) (Q 2 + 3s)s" + 7Q 2 s' 3 + 3s' 



Q 4 + (s-s') 2 + 2Q 2 (s + s') 



(12) 



n, = 



-{^ 6 [23Q 2 + 18 (-s + s')]| 

72 + 



7T 



Q 4 + (s-s') 2 + 2Q 2 (s + s') 



72 



(23 Q 2 - 18 s) (Q 2 + s) 5 + (Q 2 + s) 3 (133 Q 4 - 169 Q 2 s + 108 s 2 ) s' 
+ 2 (160 Q 8 + 6 Q e s + 3 Q 4 s 2 + 40 Q 2 s 3 - 117 s 4 ) s' 2 
+ 2 (205 Q 6 - 61 Q 4 s - 122 Q 2 s 2 + 108 s 3 ) s' 3 
+ (295 Q 4 - 37 Q 2 s- 54s 2 ) s' 4 + (ll3Q 2 -r~ 



36 s) s' 5 + 18 s' 6 



(13) 



Equation (11) corresponds to the terms containing ft-f^fi explicitly, and Eqs.(12)-(13) to the implicit 
case. The spectral function (10) contains additional terms proportional to other (independent) Lorentz 
structures, which are not written above. Collecting all three terms in (10) leads to 



ImrP(s,s',Q 2 ) 



323Q 2 + 378 (s - s') 



+ 
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4608 tt 2 4608 TT 2 [Q 4 + (s - s') 2 + 2 Q 2 (s + s')\ = 

-323 Q 12 - Q 10 (1993 s + 1237 s') - 10 Q 8 (512 s 2 + 323 s s' + 134 s' 2 ) 
+ Q 6 (-7010 s 3 + 1188 ss' 2 + 550 s' 3 ) 

+ g 4 (-5395s 4 + 7010s 3 s' + 2610s 2 s' 2 + 3146ss' 3 + 2165s' 4 ) 

- Q 2 (s - s') 2 (2213 s 3 - 2859 s 2 s' - 3099 s s' 2 - 1567 s' 3 ) 

- 378(s-s') 4 (s 2 -2ss'-s' 2 )l 



(14) 



The next step is to invoke (global) quark-hadron duality, according to which the area under the hadronic 
spectral function equals the area under the corresponding QCD spectral function. The integrals in the 
complex energy plane may involve any analytic integration kernel; this leads to different kinds of QCD 
sum rules, e.g. Laplace (negative exponential kernel), Finite Energy Sum Rules (FESR) (power kernel), 
etc. We choose the latter, as they have the advantage of being organized according to dimensionality 
(to leading order in gluonic corrections to the vacuum condensates). In this case the FESR of leading 
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dimensionality is 



/■so /■so-s /■so /■so— s 

/ ds / ds'Imn(s, S ',Q 2 ) | Hj4Z5 = / ds / ds' ImII(s,s',Q 2 ) | QCI5 . 
Jo Jo Jo Jo 



(15) 



The integration region, shown in Fig. 2, has been chosen as a triangle; the main contribution being that 

of region I, and the area included from regions II and III tends to compensate the excluded regions. Other 

choices, e.g. rectangular regions, lead to similar final results, as discussed in [6]- [7]. After performing the 

integrations, one finally obtains 

2 _ 2 s (96 Q 6 + 297 Q 4 s + 158 Q 2 s 2 - 112 s 3 ) 
l( ' ~ 9216 7T 4 (Q 2 + 2s ) A w 2 

3 ln ( Q$2^) (Q 2 + 2 5 ") ( 32 Q 6 + 67 Q 4 g + 7 Q 2 s 2 - 21 s 3 ) 
+ 9216 7T 4 (Q 2 + 2s ) X N 2 ' 

where one can recognize the standard logarithmic singularity arising from the chiral limit. In order to 
obtain the asymptotic behaviour of Fi(Q 2 ) it is essential to expand this logarithm. In fact, there is an 
exact cancellation between several terms in Eq. (16) such that the leading asymptotic term is 

^'.(O-aa^r. <") 

Qualitatively, this asymptotic behaviour agrees with expectations. 



There are two leading power corrections with no gluon exchange in the OPE of the correlator Eq.(l). The 
one proportional to the quark condensate does not contribute to Fi(q 2 ), while the other, proportional to 
the four-quark condensate, leads to 

w(p^p'\Q 2 ) = I^(± + ^)M +■■■, (is) 

where (uu) = (dd) has been assumed. The double discontinuity of this term in the (s,s') complex plane 
vanishes, so that it does not contribute to Eq. (14). 



We now turn to the extraction of F 2 (q 2 ), and consider the leading order non-perturbative power correction 
to the OPE, in this case given by the quark condensate. It turns out that the contribution involving the 
up-quark condensate vanishes (on account of vanishing traces), leaving only the piece proportional to (dd). 
The three-point function (1) becomes (see Fig. 3) 



i% 9) m (pV 2 ,Q 2 ) 

- / 



(dd) 



I 



A% ^[^{$-i)iA$-f)i a 



(k-q) 2 (k-p') 2 k 2 



■l a l u 



k 2 (k — p') 2 q 2 



(7 Q y#y ; 



') + J d 4 k 



4 ,Tr[^(^^) 7a ] 



k 2 (k — p) 2 q 2 



(19) 
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Our choice of Lorentz structure in this case is {fy M , which appears in Eq.(19), as well as in Eq.(8) where it 
multiplies i<2(<Z 2 ), but not F\{q 2 ). In fact, after some algebra 



ImII ( ^(s,s',Q 2 ) 



(dd) 



QCD 



Q 2 s' (Q 2 +3s + s>) 



1 (s'-s) 



and 



ImH"(s,s',Q 2 



Q 4 + {s- s'f + 2Q 2 ( a + s /)l 1 ^ ° 2 



HAD fe +Mw! ' r ' 



(20) 



(21) 



After substituting the above two spectral functions in the FESR Eq. (15), and performing the integrations 
one obtains 



(dd) 



% p IVI N 7I~ AJj 



2s (Q 2 + s ) + Q 2 (Q 2 + 2s ) ln( 



Q 2 



Q 2 + 2s 



(22) 



After expanding the logarithm there are exact cancellations between various terms above, leaving the 
asymptotic behaviour 

{< ' (,) -^-7TT- ) 



Q l™oo F ^ Q2} 18 Kp M N TT 2 X 2 N y Q 2 Q 4 

Qualitatively, this asymptotic behaviour does not agree with expectations. In fact, one expects F2(Q 2 ) to 
fall faster than Fi(Q 2 ) at least by a factor of l/Q [8]. Quantitatively, there is also a disagreement with 
data even at intermediate values of Q 2 , as discussed below. 



The results for the form factors Fi^(q 2 ), Eqs. (16) and (22), involve the free parameters \n and s . 
From QCD sum rules for two-point functions involving the nucleon current (2) it has been found [3]- [5] 
that Ajv — (1 - 3) x 1CT 2 GeV 3 , and ^/s^ ~ (1.1 - 1.5) GeV. The higher values of Xn and s come 
from Laplace sum rules [4], and the lower values are from a FESR analysis [5] which yields the relation 
Sq = 1927r 4 A^. After fitting Eq.(16) to the experimental data, as corrected in [9], we find An = 0.011 GeV 3 , 
and so = 1-2 GeV 2 , in line with the values discussed above. Numerically, so is well below the Roper res- 
onance peak, thus justifying the model used for the hadronic spectral function, Eq.(8). The predicted 
form factor Fi(q 2 ) is shown in Fig. 4 (solid line) together with the data, the agreement being quite good. 
A comparison of F 2 (q 2 ) from Eq.(22) with data shows a disagreement at the level of a factor two. This 
cannot be improved by attempting changes in the values of the free parameters Ajv and s , and is basically 
a consequence of the soft q 2 - dependence of F 2 (q 2 ), as evidenced by Eq.(23). 



Considering now the neutron form factors, one needs to make the change u <-> d in Eq.(2). The perturbative 
QCD spectral function, Eq.(10), involves now the combination (O3 — f2 2 ). After using the FESR Eq.(15) 
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it turns out that F\(Q 2 ) for the neutron is numerically very small and consistent with zero, except near 
Q 2 = where it diverges in the chiral limit. The explicit expression is 

F ln (Q 2 ) = — — — ^— -— ^[2 So (-75O 6 - 207Q 4 .s o -106Q 2 ,s 2 + 32 So 3 ) 

- 3 1n( * ) (Q 2 + 2s ) (25Q 6 + 44Q 4 So + 8Q 2 So 2 -6go 3 )] ■ (24) 

-\- Z Sq 

This smallness of the neutron Dirac form factor provides a nice self-consistency check of the method. Using 
F\ n {Q 2 ) ~ 0, the Sachs form factors are then proportional to F 2n (Q 2 ), which is given by 

F 2n (Q 2 ) = ^^^)[2s (Q 2 + s ) + Q 2 (Q 2 + 2 S0 ) Hq^)]- (25) 

In Fig. 5 we show the result for the electric Sachs form factor of the neutron, together with data at low 
Q 2 [10]. At higher momentum transfers, there will be a serious disagreement with experiment on account 
of the soft 1/Q 2 behaviour of F2 n (Q 2 ), Eq.(25). Since Gm(Q 2 ) for the neutron appears well fitted by 
the dipole formula, our QCD sum rule results do not agree with the data. This disagreement, though, is 
within a factor of two, i.e. not different from other recent QCD sum rule results [2]. 



In summary, Finite Energy QCD sum rules of leading dimensionality in the OPE lead to Dirac form factors 
in very good agreement with experiment for both the proton and the neutron. However, this is not the case 
for the Pauli form factor, which exhibits a soft Q 2 dependence proportional to the quark condensate. This 
is a welcome feature in several mesonic form factors where the quark condensate contributes with a 1/Q 2 
behaviour, as expected from experiment. Unfortunately, this is not the case for the nucleon. While the 
results for F 2 (Q 2 ) are dissapointing, they are not worse than those from other QCD sum rule approaches. 
In fact, the disagreement with data is within a factor two. The present method at least allows to identify 
clearly the source of discrepancy with experiment. 

We comment, in closing, on the next-to-leading order (NLO) contributions to the three-point function, 
Eq.(l), which were not considered here. On the perturbative sector we expect the gluonic corrections to 
be small, on account of the extra loop involved, plus the overall factor of a s . The NLO power correction in 
the Operator Product Expansion involves the gluon condensate. This contribution is also expected to be 
small, as it contains one more loop with respect to the leading quark condensate term. In addition, further 
suppression of about one order of magnitude would arise from numerical factors involved in the contraction 
of the gluon field tensors. On the hadronic sector, the standard single-particle pole plus continuum model 
adopted for the spectral function is well justified a posteriori from the resulting value of the continuum 
threshold s , well below the Roper resonance. 
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Figure Captions 

Figure 1. The three-point function, Eq. (1), to leading order in perturbative QCD. 

Figure 2. Triangular and rectangular integration regions of the Finite Energy Sum Rules, Eq. (15). 

Figure 3. Non- vanishing terms proportional to the down-quark condensate, Eq. (19). 

Figure 4. Corrected experimental data on Fi(Q 2 ) for the proton, [9], together with the theoretical result 

from Eq.(16) (solid line). 

Figure 5. Experimental data on Ge(Q 2 ) for the neutron [10], together with the theoretical results from 
Eqs.(24)-(25). 
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